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Abstract 

The Q 2 evolution of polarised parton distributions at small x is studied. Various 
analytic approximations are critically discussed. We compare the full evolution 
with that obtained from the leading-pole approximation to the splitting functions, 
and show that the validity of this approximation depends critically on the x —> 0 
behaviour of the starting distributions. A new analytic solution which is valid at 
small x is obtained, and its domain of applicability is discussed. 


1 Introduction 


The first moment of the polarised structure function gi(x,Q 2 ), the Ellis-Jaffe sum rule 
[1J, determines the overall spin content of the nucleon. Measurements of the Ellis-Jaffe 
sum rule [||, |3|] involve an extrapolation of gi for x —> 0, which is usually performed by 
fitting a Regge-motivated form 

g 1 (x)=Cx a (1) 

to the experimental data points with the lowest x-values. This procedure can be problem¬ 
atic, since these data points are usually taken at a relatively low scale of Q 2 ~ 1 GeV 2 , 
whereas the overall sum rule is evaluated at the average Q 2 of the experiment, which is 
typically between 3 and 10 GeV 2 . 

In the recent past, various authors have attempted to calculate the asymptotic be¬ 
haviour of gi(x) (see for example Ref. [Q for a review of the various approaches). At 
scales of low momentum transfer ( Q 2 ss 1 GeV 2 ), a non-perturbative calculation |5[ of 
the flavour singlet contribution to g\ shows good agreement with g\ at small x, but it 
should be noted that the normalisation of this non-perturbative contribution is highly 
sensitive to the only approximately known value of the vacuum quark condensate. The 
experimental discrepancy between g\ and gf in the small-x region seems to contradict 
the above result. As the singlet distribution is identical for both targets, this discrepancy 
indicates a sizeable valence-quark contribution in this region. 

With increasing Q 2 , perturbative corrections become more and more important. These 
corrections affect both the valence and the singlet contributions (A£ = J2 q (Aq + A q)) to 
gi and give rise to an evolution of the corresponding parton densities 0 
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without determining the densities themselves. These enter the above equation in the 
form of the initial distributions Aq vai (x, Qq), AE(x,Qq) and A G(x,Q%), which form the 
boundary conditions for the solution. 

In experimental measurements, these perturbative corrections are incorporated by 
rescaling the value of g\ to the average Q 2 of the experiment. This rescaling procedure 
relies on the assumption that the asymmetry gi(x) / Fi(x) satisfies exact Bjorken scaling, 
i.e. that the Q 2 -dependence of g\ coincides with that of F\. Although this assumption is 
consistent with the present data (which cover only a small range of Q 2 values at fixed x), 
there is no theoretical justification for it. In particular, examination of the polarised and 
unpolarised splitting functions || shows that g\{x)/F\{x) should indeed show only a very 
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weak Q 1 2 dependence in the large-x region, where both structure functions are dominated 
by the valence quark content, as A P gq (x) and P qq (x) are identical. In contrast to this, 
the splitting functions in the singlet sector, which dominates the small-x behaviour of F \, 
are different. The unpolarised P gq (x) and P gg (x ) have a soft gluon singularity at x — 0, 
which causes the steep rise of F\ in the small -x region. As this singularity is absent in the 
polarised splitting functions (soft gluon emission does not change the spin of the parent 
parton), one would expect the ratio j g\ (x) /F\ (x) | to decrease with increasing Q 2 . 

With the exact splitting functions 
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it is not possible to find an analytic solution to (^) with realistic boundary conditions for 
the whole range of x. By restricting themselves to small values of x (although it is not a 
priori clear which values of x can be regarded as small), various authors have attempted 
to determine the asymptotic behaviour of g\ in the limit Q 2 —> oo. One possible approach 
|| is to assume that all the Q 2 dependence is dominated by the evolution of the gluon, 
i.e. by A P gg {x). This method gives successful predictions for the unpolarised structure 
functions, due to the l/x pole in the unpolarised P gg . As this pole is not present in A P gg , 
the validity of this approach needs to be examined more carefully. 

Another possible approach flTj to the asymptotic small -x behaviour is to transform (|3|) 
into moment space and to expand around the rightmost singularity at N = 0: 

(A P) N = ^- + B + O(N) => A P(x) « A + 55(1 - a;). (4) 


This procedure yields the following approximate splitting functions^]: 
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With these simplified splitting functions, one can analytically solve (Q) for asymptotic 
values of Q 2 with realistic boundary conditions in the small-x region. This approach is 
based on the fact that the behaviour of the parton distributions at small x is governed 
by the region around N = 0 in moment space. This property can be understood from 
the ^singularity structure of the initial distributions: a logarithmic (~ 1/x ) singularity 
coincides with a pole at N = 0 in the moment transform, a power-like singularity of the 
form x a transforms into T(a + N ), which has a singularity at N = —a (see Fig. 1). ft is 
important to notice, however, that the expansion around the N = 0 pole in moment space 
agrees with the full splitting function only within a circle of unit radius (Fig. 1). Outside 
this circle, the series might still be convergent, but its value will be different from that 
given by the full splitting function. This especially affects the reliability of this approach 
for low values of a. In the extreme case a could approach —1 giving rise to a pole close 
to the boundary of the circle of convergence. 

In this letter, we examine the validity of analytical approaches to the small-x behaviour 
of g\. Section 2 contains a study of the evolution matrix on the right-hand side of (j^). Its 
properties in the case of power-like (x a ) boundary conditions are discussed, using the full 
and the leading-pole expanded splitting functions. By examining the sensitivity of the 
evolution matrix to the form of the parton distributions in the large-x region, we are able 
to assess when x can be regarded as small. In Section 3 we present an analytic solution of 
(0), which becomes exact for x —> 0. Finally, Section 4 contains some phenomenological 
implications and conclusions. 

2 Study of the evolution matrix 

Several qualitative features of the polarised parton densities can already be determined by 
inserting simple test distributions in the right-hand side of Eq. The resulting elements 
of the evolution matrix determine the local change of the parton densities with increasing 
ln(Q 2 ). Furthermore, for Q 2 jQ\ not too large one can approximate the solution of the 
Altarelli-Parisi equations by 

ag)m 2 ) = ( ag ) (l ’ G") 
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A realistic choice of test distribution is 

t(x)=x a (l-xf with (-1 < a < 0, (3 > 0), (7) 

which is similar to the analytic forms of the parton densities at used in recent fits to 
the polarised structure function data 00- The exponent a determines the behaviour of 
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the distribution in the small-a; regime, whereas the large-a; behaviour is controlled by (3. 
Variations of (3 should therefore not affect any predictions of the small-x behaviour of the 
parton distributions. This property can be used to define the range of validity of these 
predictions, i.e. to indicate if x can be regarded as small or not. 

The elements of the evolution matrix 

A ij= AP y (y)i(T) (8) 

J* y \yj 


can be computed analytically. The necessary integrals are given in Appendix A. Using 
the full splitting functions ©, we findQ 
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while the leading-pole expanded [0 splitting functions of (||) yield 
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A closer inspection of the A tJ shows that all of them diverge like x a as x —> 0. The 
behaviour in the limit x —> 0 can therefore be written as 


lim AiAx) = anx a . (11) 

1—>0 J 

Provided that both the initial quark singlet and the initial gluon distributions have power¬ 
like boundary conditions in the limit x —> 0, these most singular terms will dominate the 
right-hand side of (|2|). The replacement of the A^ 1 by the above expressions (0) in 
(@) should therefore enable us to find an analytic solution for AS(x, Q 2 ) and A G(x, Q 2 ), 
which becomes exact for x —» 0. This exercise will be performed in the following section. 

2 For simplicity, we take n/ = 3 throughout this study 
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The dij coefficients for the full and the leading-pole expanded splitting functions are 
not identical: 
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Here is the usual psi (digamma) function [|T 0 
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Figure 2 shows the A\y for a = —0.25,—0.6 and (3 = 4,9, together with the ap¬ 
proximate forms Afj and the limits a[px a . This figure displays the following important 
features of the evolution matrix in the small-x region: 


(i) Although the test distributions x“(l — x ) 4 and x a (l — x) 9 differ by less than 5% for 

.(f). 

x < 0 . 01 , the corresponding A\t differ by up to a factor of 2 in the same range. 
This clearly demonstrates that even at x — 0.01 and below the evolution is sensitive 
to the behaviour of the parton distributions in the large-;r region. The insensitivity 
of the Ay. to variations of 3 can furthermore be used to define whether x can be 

L J 1 

( f) 

regarded as small. For example, by requiring Axf- to vary by less than 30% for all 
combinations in i and j and both values of a, we fold that only x < 0.001 can 
be regarded as small, and the more conservative bound of less than 10 % deviation 
yields x < 0.0001. It should therefore be clear that the mere knowledge of g\ at 
the lowest x values accessible with fixed-target experiments is insufficient to predict 
the asymptotic behaviour of gi in the small-x limit, as the behaviour of the parton 
distributions at these values of x is still closely correlated with the distributions in 
the large-x region. 
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The convergence of the A\3 towards x a improves for smaller values of a. This 

( f) 

behaviour just reflects the fact that A\ J - contains, in addition to this leading term, 
less singular terms proportional to ln(:r). In general, these lower | A\t \. If t{x) is 
less singular that x~ l / e , the logarithmic terms are larger than the power-like terms 
for 
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where a;(a) is the branch of Lambert’s cc-function which satisfies u/(—1/e) = — 1. 
As x 0 decreases very quickly with a (x 0 ~ 10 -15 for a = —0.1), the replacement 
Alj(x) —> ciijX a , although formally still correct, loses its meaning for values of a 
close to 0 in any physically relevant region. 

(iii) While the Afj resemble the A\p for values of a close to 0, they disagree for smaller 
a. This feature becomes most striking for the A qg (see Fig. 2). The full splitting 
functions 0 predict that a positive gluon polarisation in the small-x region will 
always generate a negative contribution to the sea polarisation. In contrast, the 
leading-pole expanded splitting functions of @ predict a positive sea polarisation, 
if the gluon polarisation A G(x) is more singular than x -0 ' 5 . This behaviour can be 
inferred from the a dependence of the a %g displayed in Fig. 3. The good agreement 
for higher values of a is due to the fact that all leading contributions in ln(x) 
are contained in the N = 0 pole and hence are well approximated by the A[ l j . 
As elaborated above, these contributions remain important for a finite range in 
x > Xo > 0. The asymptotic predictions of 0 will therefore still approximate the 
full evolution, provided they are restricted to this finite range. 

(iv) The magnitude of A gg is larger by a factor 3 than the magnitude of all the other 
terms, but A gg is not more singular than any other contribution. Therefore, the 
small-x estimate of Ref. 0] is quantitative at best, and should be expected to yield 
a less accurate prediction than the corresponding estimate of the unpolarised dis¬ 
tributions. 
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Table 1: Residues of the polarised splitting functions in A-mo merit space. The residues 
for all negative integers with N < —2 are identical. 


(v) The agreement between leading pole expanded and full splitting functions is better 
for the A gq and A gg than it is for A qq and A qg . This feature can be understood 
from the relative magnitude of the residues in the corresponding splitting functions 
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(Table [l]): the N = 0 residue is dominant only in the P gq and P gg splitting functions, 
the other two splitting functions contain residues for N < 0, which are twice as big 
as the N — 0 residue. 

It should be clear from the above that the leading-pole expansion of Ref. 0 gives a 
reliable approximation to the evolution matrix in the small-X region, provided that the 
initial distributions are significantly less singular than x l ' e . For more singular distribu¬ 
tions, this approach results in a manifestly different evolution matrix and hence will yield 
a different small-x behaviour of the polarised parton distributions. 


3 Solution of the Altarelli-Parisi equations in the limit 

x —> 0 

Provided both polarised singlet quark and gluon densities have power-like boundary con¬ 
ditions in the small-x region, 

AE(x, Qq) ~ x aq , AG(x, Ql) ~ x aa with — l<a 9 ,a G <0, (14) 

one can find a solution of the Altarelli-Parisi equations which becomes exact in the limit 
x —> 0 and has the form 


A q va i(x,Q 2 ) = R v (Q 2 ,Ql) x av 

A E(x,Q 2 ) = R qq (Q 2 ,Ql)x aq + R qg (Q 2 ,Ql)x aa , 

A G(x : Q 2 ) = R gq (Q 2 ,Ql)x^+R gg (Q 2 ,Q 2 0 )x aG . (15) 

A detailed derivation of this solution and the explicit forms of the /?-functions is given in 
Appendix B. 

The above bounds on a cover the whole theoretically allowed range: as the first 
moments of the distributions have to be finite, we find a > — 1. Furthermore, inspection 
of the singularity structure of the evolution equations (Fig. 1) shows that any initial 
distribution, which is finite in the small -x region, will develop a logarithmic divergence 
due to the N = 0 singularity of the splitting functions. The case of finite or logarithmic 
boundary conditions can be treated correctly with the leading-pole approximation - its 
asymptotics are discussed in [0. In a previous analysis 0 of the experimental data on 
polarised structure functions we have found a q = a v ~ —0.55. The experimental data 
used in this analysis were insufficient to determine ct G , and therefore it was fixed to be 0. 
In contrast to this, more recent measurements at lower Q 2 || favour ct G < 0. 

As we have neglected all contributions of order ln(x) in the above solution, we expect 
it to be reliable only for x < Xo(max(a g , a G )). In order to compare this approach with the 
leading pole expansion of Ref. [0 and the numerical solution of (||) with the full splitting 
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functions, we have evaluated the distributions for Ql = 4 GeV 2 and Q 2 = 100 GeV 2 , 
using rif — 3,A® cd = 200 MeV and the following initial distributions: 

AS(i, Qq) = N q x aq {l — x) 13 

A G(x,Ql) = N G x aG ( l-xf (16) 

A q val (x,Ql) = N val x av ( l-x) 13 . 

To illustrate the validity of the various approximations, we adopt the following parameter 
values: a q , a G , ot v = —0.6, —0.25, (3 = 4, 9, and for simplicity we take N q = N g = N v = 1. 

Figures 4 (a), (b) and (c) show the behaviour of the gluon, singlet quark and va¬ 
lence quark distributions respectively, at small x and Q 2 = 4, 100 GeV 2 . The initial 
distributions x a {l — x) 4 are indicated as solid lines. 

Starting with the gluon distribution (Fig. 4(a)), we see that for x < 10~ 2 , the leading- 
pole approximation to the splitting funcitons (dotted lines) gives excellent agreement with 
the full evolution (dashed line), especially for values of a q , a G close to 0. This is consistent 
with the agreement between the corresponding A gg functions shown in Fig. 2 and can be 
understood from the N — 0 dominance in the A P gq and A P gg splitting functions. In con¬ 
trast, the x a approximation (short-dashed line) significantly overestimates the evolution 
in the x range shown, espcially for a q , a G close to 0. Convergence of this approach can 
only be observed at even smaller values of x. Note, however, the sensitivity to the large-x 
behaviour. While both the dotted and the dashed lines are computed with (3 = 4, the 
dot-dashed curve corresponds to full splitting function evolution for (3 — 9, i.e. a softer 
large-x distribution. Evidently there is a significant sensitivity to the behaviour at large 
x even for x values as small as 0(1CT 3 ). This casts doubt on the idea of using data on the 
evolution of the small-x structure functions alone to determine the gluon distribution. 

For the singlet quark distribution (Fig. 4(b)) the situation is rather different. Here 
the leading-pole approximation overestimates the evolution at small x. This is readily 
understood from the behaviour of the corresponding A qq and A qg functions in Fig. 2, both 
of which are systematically more positive for the leading-pole splitting functions. In fact 
we see that for a q = —0.25 and a G = —0.6, the full evolution gives a negative singlet 
distribution at small x, whereas the leading pole splitting functions give a positive distri¬ 
bution. Notice also that the evolution is less sensitive to the large-x behaviour (compare 
the dashed and dot-dashed curves which correspond to (3 — 4, 9 respectively) than for the 
gluon distribution. For a q = a G = —0.6, the x a approximation is quite reasonable, and 
certainly better than the leading-pole approximation. However the opposite is true when 
both a q , a G are close to 0. 

Finally, Fig. 4(c) compares the valence quark evolution in the various approximations. 
This depends only on A P qq , and so the behaviour here is a direct reflection of the corre¬ 
sponding A qq shown in Fig. 2. In particular, for a q = —0.6 the x a approximation is very 
good, while the leading-pole approximation overestimates the evolution at all x values 



shown. For less singular small-x behaviour (a q = —0.25), however, both approxima¬ 
tions reproduce the full evolution, the leading-pole approximation showing slightly better 
convergence for x > O(10 -4 ). 

In practice, the normalisations of the singlet quark and gluon distribution, N q and Nq, 
will not be the same. As the evolution of the gluon densitity is dominated by the gluon- 
to-gluon splitting, it will be almost unaffected by changes of N q . Only if N q is one or more 
orders of magnitude larger than Nq, will the impact of quark-to-gluon splitting become 
visible. More drastic effects of a change in the relative normalisation can be expected for 
the quark singlet distribution, as contributions from quark-to-quark and gluon-to-quark 
splitting have the same magnitude but opposite signs (cf. Fig. 3). Therefore, a relative 
increase of Nq yields a faster evolution of the quark distribution to negative values. 

The convergence properties of the different analytic approaches are almost unaffected 
by changes in the normalisation. Only for Nq N q do we find that convergence of the 
x a approximation to the singlet distribution sets in for smaller values of x. This simply 
reflects an increased impact of the gluon-to-gluon splitting. 


4 Conclusions 

In this letter we have studied the feasibility of two different analytic approaches to the 
evolution of polarised parton densities at small x, finding that none of these approaches 
is able to give reliable predictions for the whole theoretically allowed range of boundary 
conditions in the small-x region. In the leading-pole expansion [|7j. |[], the full splitting 
functions A P iq are replaced by the leading terms of their Laurent series around N = 
0 . As this approach correctly reproduces all terms proprotional to In x generated in 
the evolution, it is found to be in good agreement with the full evolution if the initial 
quark and gluon distributions are less singular than x~ l ! e . For more singular boundary 
conditions, only the gluon distribution is reproduced correctly, in particular the quark 
distribution is overestimated. Keeping only terms with powerlike singularities in the 
evolution equation, we were able to derive an exact solution of this equation in the limit 
x —> 0. As we have neglected all logarithmic terms in this approach, its convergence is best 
for boundary conditions of quark and gluon distributions more singular than x _1 / e . For 
less singular boundary conditions, this approach still converges towards the full solution, 
but its predictions are far away from the full solution for any realistic experimental value 
of x. 

We have also shown that the evolution of the polarised gluon distribution is sensitive 
to the shape of this distribution in the large-x region. This observation raises doubts 
on the possibility of determining the gluon polarisation from the evolution of g\ in the 
small-.T region. It furthermore demonstrates the need for complementary measurements 
of A G(x) (e.g. from J/T-production or direct-y measurements). 
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We have seen that the effects of the evolution on the quark distributions in the small-x 
region are rather small, as the quark-to-quark and the gluon-to-quark splitting contribute 
with opposite signs. The gluon distribution is indeed rising with increasing Q 2 , but 
only contributes to g\ at order a s (Q 2 ). Bearing in mind that AG contributes with a 
negative coefficient function to g±, one expects that g\ will become negative at small x for 
asymptotic values of Q 2 , due to the gluonic contribution and the negative sea polarisation 
generated from g —> qq splitting. 

In general, the effects of the evolution on the polarised parton densities will be more 
moderate than the effects on the unpolarised densities. The assumption of approximate 
scaling for gi (x)/F\ (x) in the small -x region is therefore rather doubtful. It seems more 
realistic to assume approximate scaling for g±(x) for the range of fixed-target experiments, 
due to the partial cancellation of quark and gluon evolution as explained above. 


Acknowledgements 

Financial support from the UK PPARC (WJS), and from the Gottlieb Daimler- und Karl 
Benz-Stiftung and the Studienstiftung des deutschen Volkes (TG) is gratefully acknowl¬ 
edged. This work was supported in part by the EU Programme “Human Capital and 
Mobility”, Network “Physics at High Energy Colliders”, contract CHRX-CT93-0357 (DG 
12 COMA). 


A Convolution integrals of the test distribution 

For the test distribution 


t(x) = x a {l — x) 13 with (—1 < a < 0, /3 > 0), (17) 

the convolution integrals 

Aj = t — APiM t (-) ( 18 ) 

A y \yj 

on the right-hand side of (||) can be expressed in an analytic form. From the explicit 
forms of the splitting functions given in (||) and (Q), one sees that the required integrals 
are 


Ai(x) 


r 1 dy 1 (x \ a ( x 
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All these functions diverge like x a as x —> 0, and the leading singular behaviour at small 
x is found to be 


Ai(x) ^-4 x a [-^(-a) - 7 e] 
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A 4 (x) x a . 
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B Analytic solution of the Altarelli-Parisi equations 

for x —> 0 

A solution of the Altarelli-Parisi evolution equations (|2|) can never be more singular at 
x = 0 than the starting distributions. It follows that the most singular parts of the valence 
quark, singlet quark and gluon distributions can be obtained by inserting the following 
ansatz 

Aq va i(x,Q 2 ) = R V (Q 2 , Ql) x av 
A £(x,Q 2 ) = R qq (Q 2 ,Q 2 0 )x^+R qg (Q 2 ,Q 2 0 )x aG , 

AG(x,Q 2 ) = R gq (Q 2 ,Ql)x aq + R gg (Q 2 ,Q 2 0 )x aG (21) 

into (2j). Keeping only terms proportional to x a ", x aq and x a ° on the right-hand side, we 
obtain the following evolution equations for the R coefficients {(3q = 11 — 2/3 rif) : 

-^—R v (Q\Ql) = -l aqq (a v )R v (Q 2 ,Q 2 0 ) 
a In a s p 0 
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As we are interested in the asymptotic solution for the full splitting functions, all a,-, in 
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the above are a); . 
Introducing 

s = In 
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/ln(g 2 /A 2 ) \ 


\ln(Ql/A 2 )) 
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the general solution of these equations reads 
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Po 


Po 


where the Rij±(Qo) are determined by the boundary conditions at Qf. As we assume that 
the initial distributions for the quark singlet and the gluon have the form 

AS(x, Ql) = N q x a \ A G(x, Q 2 0 ) = N G x aG , 


(25) 


these constants are determined to be 
UJ +( a q) a ggi a q) 


n (n 2 \ — UJ +( a q) a ggi a q ) * r d /p2\ _ U-(a q ) a qq (a q ) , r 

o;-(_(a g ) — LU-(a q ) a; +l Q; (}J 

i? OT+ (g 2 ) = a ^ a i) \ Nq, Rgq-iQl) = - a 9 q(<Xq) N 

9q y u/ u+\ot q ) -cn_(a g ) q: 9q u + (a q ) - u_(a q ) q 

D AO 2 ! — _ a, ig ( Q 'G ) _/y 73 ('O 2 ') —_ a qg{ a G) _yy 

99+WoJ cn+(a G )-cu_(a G ) 9 ’ - u+ (a G ) - u_(a G ) 9) 

73 in 2 '! _ ^+( Q 'g) ~~ a qg( a G) ]\T TD (r>2\ _ _ lJ-(Q'g) — Cl qq {c»L G ) M 

Kgg+Wo) w+ ( aG )- w _(a G ) J b> “gg-Wo) u+{a G ) - uj.(a G ) 9 ' 


( 26 ) 
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Figure Captions 

Figure 1 Singularity structure of the evolution equations in the complex TV-moment 
plane. Dots (•) denote the poles of the splitting functions, and the cross (x) indi¬ 
cates the small-x singularity of the initial distribution. The leading-pole expansion 
only converges to the splitting function in the unit circle around the origin. 

Figure 2 Elements of the splitting matrix for the test distribution x a {l—x)P. Solid line: 
full splitting functions for fi — 4, long-dashed line: same for (3 = 9, short-dashed 
line: most singular x a contribution, dotted line: leading-pole expanded splitting 
functions for f3 = 4, dot-dashed line: same for (3 — 9. For better visibility, all 
elements are multiplied by x. 

Figure 3 Coefficients of the most singular pieces in the splitting matrix for the full 
(left) and the leading-pole expanded (right) splitting functions. Solid line: a qq , 
long-dashed line: a qg , short-dashed line: a gq , dotted line: a gg . 

Figure 4 Evolution of test distributions for gluons (xAG(x, Q 2 )),(a), singlet quarks 
(xAS(x, Q 2 )),(b) and valence quarks (xAq va i(x,Q 2 )),(c) as described in the text. 
Solid line: starting distribution at 4 GeV 2 , long-dashed and dot-dashed line: evolved 
distributions at 100 GeV 2 for different large -x behaviour at Q%, short dashed line: 
result of x° approximaton, dotted line: result of leading-pole approximation. 
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